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We propose and analyze difTerent schemes to probe the quantum nature of nanoelectromechanical 
systems (NEMS) by a tunnel junction detector. Using the Keldysh technique, we are able to 
investigate the dynamics of the combined system for an arbitrary ratio of eV/h^l, where V is the 
applied bias of the tunnel junction and fl the eigenfrequency of the oscillator. In this sense, we 
go beyond the Markov approximation of previous works where these parameters were restricted 
to the regime eV/hO, 3> 1. Furthermore, we also go beyond the Born approximation because we 
calculate the finite frequency current noise of the tunnel junction up to fourth order in the tunneling 
amplitudes. 

Interestingly, we discover different ways to probe both position and momentum properties of 
NEMS. On the one hand, for a non-stationary oscillator, we find a complex finite frequency noise 
of the tunnel junction. By analyzing the real and the imaginary part of this noise separately, we 
conclude that a simple tunnel junction detector can probe both position- and momentum-based 
observables of the non-stationary oscillator. On the other hand, for a stationary oscillator, a more 
complicated setup based on an Aharonov-Bohm-loop tunnel junction detector is needed. It still 
allows us to extract position and momentum information of the oscillator. For this type of detector, 
we analyze for the first time what happens if the energy scales eV, HQ, and fcsT take arbitrary 
values with respect to each other where T is the temperature of an external heat bath. Under these 
circumstances, we show that it is possible to uniquely identify the quantum state of the oscillator 
by a finite frequency noise measurement. 

PACS numbers: 85.85. +j, 72.70.+m, 07.50. Hp, 42.50.Lc 



I. INTRODUCTION 

Nanoelectromechanical systems (NEMS) have become 
a promising playground for probing the quantum be- 
havior of mesoscopic objects, theoretically as well as 
experimentally-^. The diverse reasons to study NEMS 
are their vast number of (possible) applications as for in- 
stance the measurement of mass, force and positiorP^^ 
with high precision. 

Nanomechanical systems being at the boarder of clas- 
sical to quantum are also being studied from a very fun- 
damental point of view. This includes the observation, 
measurement and control of quantum states of a meso- 
scopic mechanical continuous variable system such as a 
harmonic oscillator. Superconducting qubits electrically 
coupled to the mechanical system have been successfully 
used to characterize the mechanical resonator's quantum 
stateP. 

Making quantum effects visible in nanomechanical sys- 
tems calls for ultralow temperatures and low dissipa- 
tion. The goal of observing the quantum mechani- 
cal ground state of a harmonic oscillator requires tem- 
peratures ksT <C hil. Recently, this goal has been 
achieved by using a microwave-frequency mechanical os- 
cillator, with a frequency of / « 6 GHz which allowed 
cooling to the ground state with conventional cryogenic 
refrigeratiorP. Further proposals of cooling a nanome- 
chanical reson ator coupled to an optical cavity have been 
proposecP^'^ and experimentally implementecP^^. 

The theoretical treatment of NEMS widely uses a 



Markovian master equation approacfP^''^ with a few 
exceptions, for instance, the work by Wabnig et al^ 
and Rastelli et al^ where a Keldysh perturbation the- 
ory has been employed. Here, we also make use of 
the Keldysh technique because it allows us to treat the 
non-equilibrium system fully quantum mechanically and, 
furthermore, to carefully investigate the non-Markovian 
regime where eV <^ Ml. Since we are interested in the 
quantum nature of the oscillator, it is important that the 
temperature T and the applied bias V of the tunnel junc- 
tion are not much larger than the eigenfrequency of the 
oscillator. Otherwise, the oscillator would be heated and 
low energy properties inaccessible. 

The article is organized as follows. Our key results 
are summarized in Sec. [llj In Sec. |III[ we introduce the 
generic model. This is followed by an introduction of the 
formalism we use in Sec. |IV| with subsections focusing on 
the fermionic reservoir and on the oscillator dynamics. 
The main part of this article is presented in Sec.jVj where 
we discuss the calculation as well as the results for the 
finite frequency current noise. Finally, we conclude in 
Sec. ED 



II. KEY RESULTS OF THE PAPER 

The motivation of our work is to study an experi- 
mentally feasible setup in which the quantum nature of 
NEMS can be probed by current noise measurements 
of a tunnel junction detector. A quantum NEMS can 
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be described by a quantum harmonic oscillator which is 
a continuous variable system characterized by two non- 
commuting operators x and p. Therefore, it is desirable 
to have a detector at hand that can measure expecta- 
tion values with respect to x-dependent observables, p- 
dependent observables, as well as observables that de- 
pend on both X and p. 

In Ref. Ifji Doiron et al. have proposed a setup which 
could be used for position and momentum detection of 
NEMS. This setup consists of two tunnel junctions form- 
ing an Aharonov-Bohm (AB) loop. There, it is possible 
to tune the relative phase between the tunnel amplitudes 
(where one depends on x and the other one not) via a 
magnetic flux penetrating the AB loop, see Fig. [T] for 
the schematic setup. In such a setup, the symmetrized 
current noise 



^symH = I ydte''^*({A/(t),A/(0)} 



(1) 



(with the current fluctuation operator A/(i) = i{t) — 
(/)) of the tunnel junction detector can either contain 
information on the oscillator's position spectrum 



S'a;(w) 



dte-* {{x{t),x{0)}), 



(2) 



then S. 
spectrum 



sym('^) 



5*2; (w), or the oscillator's momentum 



Sp{uj) 



dte^-' ({p(i),p(0)}). 



(3) 



sym (^) 



Sp{uj), see also Eq. (55) below. The 



then St 

current noise Ssym{^) can also contain information on 
both, the position and the momentum of the oscillator. 
The former case has been coined cc-detector and the latter 
case p-detector. We note that S'sym(w), Sx{<^), and Sp{uj) 
are properly defined above in Eqs. ([l])-(3| for a stationary 
problem. In the non-stationary case, which is also subject 
of discussion in this work, these quantities do not only 
depend on a single frequency lo but on one frequency 
argument and one time argument instead, see Eq. (39). 

To be more specific, in this article, we call an x- 
detector, a detector that allows to measure expectation 
values of the oscillator's position operator x, i.e. (x), 
{xx), etc.. Similarly, we call ap-detector, a detector that 
allows to measure expectation values of p, the oscilla- 
tor's momentum operator, i.e. {p), (pp), etc.. In Ref. IT6l 
switching from the x-detector to the p-detector is then 
accomplished by tuning the relative phase between the 
tunnel amplitudes. The main difficulty of this setup is 
the need of long coherence times and length in the AB 
loop to make the switching possible. Here, we show that 
the AB setup can be avoided. We find that the current 
noise of the coupled oscillator-junction system with one 
tunnel junction only, already can be used for momentum 
detection due to the complex nature of the current noise 
when the oscillator is in a non- stationary state. This is 



the first key result of our work, specified and intensively 
discussed in Sec. IV B 2l below. 

We further investigate the current noise stemming from 
a stationary oscillator up to fourth order in the tunnel- 
ing amplitudes, thereby going beyond the Born approx- 
imation. Most importantly, we extend previous results 
of Ref. to the non-Markovian regime without any re- 
strictions on the relative magnitude of the energy scales 
eV, hfl, and fc^T. We show that peaks in the finite 
frequency current noise at uj — ±il (both for the x- 
detector and the p-detector) are a fourth order effect. 
In the Markovian regime, the peaks in the position de- 
tector signal are always much larger than the ones in 
the momentum detector signal. This is different in the 
non-Markovian regime. There, we even find a larger sig- 
nal for the momentum detector compared to the position 
detector, clearly demonstrating that the non-Markovian 
regime is the preferred regime to operate the momentum 
detector. The detailed understanding of the x-detector 
and the p-detector developed in this article allows us to 
uniquely identify the quantum state of the oscillator by 
a finite frequency noise measurement. This is the second 
key result of our work, specified and intensively discussed 
in Sec. lVC4l and Sec. ITCSl below. 



x{t) 





FIG. 1. Schematic setup for the realization of a position 
detector which can be extended by the parts enclosed in the 
red dashed box to a momentum detector. The total tunnel 
amplitude in the case of the position detector is given as /3 = 
P{x) — to + tix and in the case of the momentum detector as 
/? = P{x) + I3 = to + tie^^x, where the relative phase rj between 
to and ti can be tuned via a magnetic flux $ penetrating the 
AB loop. If the oscillator is in a non-stationary state, already 
the parts without the elements in the red dashed box serve as 
a position as well as a momentum detector. 



III. MODEL 

The system we consider consists of a nanomechanical 
harmonic oscillator coupled to a biased tunnel junction. 
In Ref. [19] an experimental realization is shown, where 
electrons can tunnel from an atomic point contact (AFC) 
onto a conducting oscillator. The coupled system is de- 
scribed by the following Hamiltonian 



(4) 
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with 



and similarly 
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(5) 
(6) 



(9) 



where ifosc describes the oscillator with x and p being the 
position and momentum operator of the oscillator with 
mass TO and frequency f7, respectively. H^^s contains the 
fermionic reservoirs of the left and right contacts. The 
oscillator couples to the tunnel junction via the tunneling 
Hamiltonian 



P c\ Cr -\- h.C. 

l,r 



(7) 



Here (cj) annihilates (creates) an electron in reservoir 
i = l,r. Motivated by the experimental setup in Ref. UHl 
we take the oscillator to act as one of the fermionic reser- 
voirs. Therefore, the tunneling gap depends on the po- 
sition of the oscillator, modifying the tunneling ampli- 
tude of the APC. For small oscillator displacements x, 
we assume linear coupling of the oscillator to the tunnel 
junction with a tunnel amplitude /3i. Hence we obtain 
/3 = [/3o + /?! x] , with I3q being the bare tunneling ampli- 
tude. Here, we allow for complex tunnel amplitudes (3q 
and /3i as previously discussed in Refs. HHI and HOI With 
ry we denote the relative phase between the tunnel am- 
plitudes, i.e. we write /Sq — Iq and /3i = tie*'' where 
to,ti G ^. A possible experimental realization of the fi- 
nite and tunable phase ij is discussed in Ref. [161 As a 
consequence this phase 77 gives rise to the possibility to 
detect the oscillator's momentum expectation value (p^), 
present in the current noise. 



IV. GREEN'S FUNCTIONS OF THE 
FERMIONIC RESERVOIR AND OF THE 
OSCILLATOR USING THE KELDYSH 
TECHNIQUE 

The true non-equilibrium, non-Markovian quantum 
behavior of the coupled system is the subject of our 
interest. Th erefore, we make use of the Keldysh 
formalisrEpUm for our calculation. The quantities being 
accessible in the experiment are for instance the tunnel 
current and the current-current correlator, the noise of 
the tunnel junction. These will also be the main objects 
of interest in this article. We employ a perturbation the- 
ory in the tunnel Hamiltonian Htun and calculate the 
noise up to fourth order in the tunneling. 

The current operator is given by J = —eNi, where 
Ni = cjci counts electrons in the left reservoir. We then 
write the current operator as 



I = e jo + xji 



(8) 



with ji 



and hi — Tt +7^ . The operator % 



is given by 7i = J2i.rf^i^l^r with i e {0, 1}. 



A. Reservoirs Green's functions 

The fermionic Green's functions of the left and right 
reservoirs (free electron gas) Gi^r, are given on the 

Keldysh contour C by Gi^r{t,t') = -i (^T^ Q,r(t)c3'^(t'))- 
Tc denotes the time ordering operator on the Keldysh 
contour, placing times lying further along the contour to 
the left. Figure [2] shows the Keldysh contour C. The con- 
tour consists of a lower branch, C_ on which time evolves 
in forward direction and of an upper branch, C+ where 
time evolves in backward direction. Switching from times 
lying on the contour to real times is done by analytic con- 
tinuation. The Keldysh Green's functions Gi^r{t,t') can 
then be represented by a matrix. A Fourier transforma- 
tion leads to the following Green's functions 



Gf+ 



(c.) 



ni^ri^) - 1/2 
nir{^^) - 1 



(10) 



27ripo 



ni^ri'-^) — 1/2 



Here we made use of time translation invariance and 
assumed a constant density of states in the left and 
right reservoir pi = pr ~ po- The applied finite bias 
fJ-r ^ fJ-i = eV is included in the Fermi distribution func- 
tions ni = n{uj - eV/2) = [exp(/3(a; - eV/2) + 1]"^ and 
Ur = n{uj + eV/2) = [exp(/3(w -I- eV/2) + 1]"\ The in- 
verse temperature of electrons in the reservoirs is /3 = 
l/fcsT and we use units where h= 1. 



c 



FIG. 2. Keldysh contour C with the lower branch C_ and 
the upper branch C+. 



B. The oscillator 

Since the oscillator modulates the tunneling of elec- 
trons and therefore has impact on the measured average 
current and current-current correlator, it is important to 
understand the significance of the oscillator's state. We 
distinguish between an oscillator in a stationary state and 
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one in a non-stationary state. We justify this differenti- 
ation by arguing that for short times after the measure- 
ment, the oscillator will certainly be non-stationary. The 
dominating timescale here is the one given by the oscilla- 
tor itself, 1/ri, which has to be compared to times scales 
on which the damping of the oscillator due to the tunnel 
junction and the external heat bath happens. In the non- 
stationary case, we cannot make use of time translation 
invariance in the oscillator's correlation function D(t, t'). 
For longer times however, the assumption of stationar- 
ity is justified since the oscillator can equilibrate with 
the environment and reach a steady state. The oscilla- 
tor's correlation function now only depends on the time 
difference t — t'. 

We work in the following with the oscillator operators 
given in the Heisenberg picture as x{t) = xcos{^t) + 
p/(mn) sin(J7t) and p{t) — pcos{iU) — x{'m,n) sm{nt). 
We also define the aforementioned oscillator correlation 
function D{t,t') in Keldysh space as 



where we defined 



D{t,t') = -i{T,x{t)x{t')) 



(11) 



When we later investigate the second order noise we 
consider both the stationary situation and the non- 
stationary situation. The following relation then is a very 
useful one 

x{t + t') = i(t') cos(m) + sin(m) . (12) 

mil 

For calculations up to second order, we look at the in- 
fluence of stationary and non-stationary oscillator states 
on the current noise, in fourth order we restrict ourselves 
to the stationary case. Hence, we are interested in a 
clear definition of the oscillator's correlation functions 
and spectral functions in the stationary case, which will 
be addressed now. 



1. Oscillator correlation functions in the stationary case 

Considering the stationary case, we give useful ex- 
pressions for the oscillator's correlation functions which 
later allow us to identify the oscillator's power spec- 
trum in X denoted by Sx{u>) and in p denoted by 5*^(0;). 
From Eq.(ll) the correlation function where t G C+ and 



t' £ C_ is given by 

iD+-{t,t') = {x{t)x{t')) 



= -1x1 cos{n{t - t')) + sininU - t')) 
2 \ m\l 



+ xpLco&{VL{t + t')) + 



mil 



sin(17(t + t')) 



(13) 



and similar for t G C_ and t' G 
iD-+{t,t') = {x{t')x{t)) = 

[P, 



= -(xl cos(fl(t - t')) - sin(n(t - t')) + 

2 \ mil 



+ x'i cos{n{t + t' j) + 



{P,x} . 



(14) 



-2 , P' 
X ± 



(15) 



and [•, •] denotes the commutator and {•, •} the anti- 
commutator. Since, here we deal with the stationary 
case, the expectation values (S?.) and {{p, x}) appearing 
as prefactors of functions depending on t + t' equal zero 
which one can easily check by using any stationary state, 
e.g. number-states. As one would expect, the correlation 
function now is a function of the time difference t — t' 
only. The Fourier transform of the correlation functions 
then yields 

^D+-{u) = l{xln+{u, n) + ^7^-(c., n)) , (i6) 
^D-+{u) ^ ^{xln+{u;,n) - ^7^-(c.,^l)) . (i?) 

Additionally, we introduce the two momentum correla- 
tion functions iP+-{t,t') = {p{t)p{t')) and iP-+{t,t') = 
{p{t')p{t)). The same arguments as for iD^^ {t,t') lead 
here to the following Fourier transforms 



iP+-{w) = ^(^plTz+{ij,n) + mni[p,x]TZ-{w,n)y 



iP-+{u) = ^{p\n+{iu,9) ^ mm[p,x]n-{u:,^l)) , 



(18) 
(19) 



where similar to above 



p^ — m il x ± p . 
We introduce the functions 7?.^(w, Vl) as 



(20) 



n+{Lo,n) = Ti[5{Lo + n) + 5{Lo-n)], (21) 



(22) 



The coupling of the oscillator to two environments, 
namely an external heat bath and the tunnel junction, 
being at the temperatures Tcnv and fc^Tjunc — eV/^r^ 
respectively, introduces a damping of the oscillator with 
damping coefficients 70 and 7+ respectively. The oscil- 
lator dynamics due to the coupling to the tunnel junc- 
tion can be calculated by solving a Dyson equation for 
the oscillator correlation function D{t,t') where the self- 
energy is taken to lowest non-vanishing order in the 

tunnel Hamiltonian, i.e. S(t,t') = iTchi{t)hi{t')^ . 

Using the Keldysh technique as done in Refs. |T7] and 
l24l the oscillator dynamics and the damping coefficient 
7+ = TT/jgi^/m can be calculated. The coupling to the 
external heat bath can be added phenomenologically, 
or particularly as an interaction with a bath of har- 
monic oscillators. The total damping then follows as 
7tot = 7o + 7+ • We can assign an effective temperature 
Teff to the oscillator with ^totTcS = 7+ Tjunc + loT^nv 
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This leads to the general case for TZ^ by replacing the 
i5-functions in Eqs. ( 21|22 ) by a Lorentzian, where we 
include both sources of damping and an oscillator fre- 
quency Q — > -^/il^ — 72. For this damped case we can 



write for 7?.+ (w, fl) and TZ (w, fi) 



47tQt uj - 7^2^^ 



(23) 
(24) 



We want to stress that we have not made any assumption 
on the initial time as e.g. t' = 0. With the above made 
definitions, the following relation (which has to hold for 
any bosonic correlation function in Keldysh space) 



(25) 



can easily be verified, here i and j are Keldysh indices. 
This concludes our discussion on the oscillator correlation 
functions. We now turn to the spectral functions Sx{ijj) 
and Sp{uj). 



We can also write down the spectra using the Keldysh 
Green's function (uj) = (lu) + D '"(w) which 
yields 

SAoj) = liD''{Lj), (29) 
Sp{uj) = ^m^n^ iD^{uj) . (30) 

To further simplify the notation, we introduce 

Q{uj)^'-[D+-{u;)-D-+{lo)] = 

= ^7^;(.,^7), (31) 

which is used later in the fourth order noise calculation. 



2. The oscillator's spectra in the stationary case 

The symmetrized power spectrum, in general defined 
as 5/ die*"* (^{f{t),f{f)}'^ of the oscillator quantities 

T = a;,p is an observable that can be measured by e.g. 
current noise measurements (as discussed below). Both 
Sx{uj) and Sp{uj) can be measured through the current 
noise S'syin(w). The expressions for these power spectra 
are given by 

5,(w) = ^ J dte^-' {{x{t).i{t')}) = 

= dte'^U [D+-{t,t')+D-+{t,t')] = 

= \ [D+-{lo) + D-+{u^)] = 

= \{x^)n+{i^M). (26) 

and 

Sp{u) = \ j dte^-'{{p{t),p{t')})^ 

= \j dte^-U [P+-{tX) + P-+{tX)\ = 
= iz [P+-(C.)+P-+(C.)] = 

= i(p2)7^+(^,^7). (27) 

The momentum and position spectrum are related via 
the relation 

Sp{ijj) = rn^Vi^ Sx{u) . (28) 



V. CURRENT NOISE CALCULATIONS 

In this section, we cover a variety of aspects when deal- 
ing with the current noise. For all different aspects we 
find expressions for the noise which are valid for an ar- 
bitrary rj and therefore include the x-detector as well as 
the p-detector. 

The first part is dedicated to the noise in second or- 
der perturbation theory, where we furthermore make 
the distinction between a stationary harmonic oscillator 
and a non-stationary one. Beside the Markovian regime 
{eV ^ Ml), the Keldysh formalism also allows us to in- 
vestigate the non-Markovian regime [eV ^ Ml). The 
main results in this section are that the current noise for 
a non-stationary oscillator can in principle be complex. 
In this detectable complex noise would allow for 

a nearly complete determination of the oscillators covari- 
ance matrix cTy = tr (/5{T), Tj}/2), where T = {x,p)'^. 
The covariance matrix allows for a complete description 
of the oscillator's quantum state. 

For the stationary oscillator we recover a noise that is 
real and in accordance with the Wiener-Khinchin theo- 
rem. This noise is the well know noise of a bare b iased 
tunnel junction which shows kinks at |a;| = mod- 
ified by the oscillator leading to kinks at \Ld\ = \ V ztn\. 

In the last part of this section, we deal with the noise 
up to fourth order in the tunneling amplitudes. Then 
we restrict ourselves on the stationary case. On the one 
hand, the fourth order contributions modify the kinks, 
on the other hand, they give rise to resonances stemming 
from the oscillator correlation functionJSEHIII, 

We now introduce the perturbation theory leading to 
the noise expression. 
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A. Overview 

In general, the expression for the current-current cor- 
relator in the Keldysh formalism is given by 

'Jirs)) (Hr,)) , (32) 



where a (f2, V) is given in Eq. (43). Our result for the 



since we consider only the second and fourth order cur- 
rent noise we write 



S{T3,Ti) = S'^^\T3,n) + S^^\T3,n) , (33) 



where 



and 



(34) 



1 



dTidT2 ( Tc -fftun(Tl)i/tun(T2)/(T3)/(r4^ 



dTidT2 (TciJtun(Tl)/(r3) ) ( Tci/tun (t2 )/(t4 



(35) 



The general expression for the average current is given 
by 



(/(t)) = (r,e 



m 



(36) 



To second order in the tunneling amplitudes, the average 
current can be calculated by 



{I{t)) ^ -I j dT {T, Htnn (T)/(t) ) , (37) 

and we obtain 



+2cos{T])totieV {x{t)) + sin{T])tot: 



(Pit)) 



+tleV{xit)m~^'^-i^,V)], (38) 



average current is in accordance with Ref. [161 The cur- 
rent noise we calculate is always the frequency-dependent 
(and in the non-stationary case also time-dependent) 
symmetrized current noise, defined as 

5sy„,(c^,i') = ^J dte'^' [S-+{t,t') + S+-{t,t')] (39) 

where for S ^{t,t'), t E C_ and t' e C+ and similar 
for 5*^ {t,t'), here C_ and C+ are the lower and upper 
branch of the Keldysh contour C, respectively, see Fig. [2] 
B. Current noise to second order in the tunneling 
amplitudes 

We now turn to the calculation of the current noise 
to second order. With the current operator / already 
being first order in the tunneling amplitudes, the current 
noise in the Born approximation is given by the following 
expression 



(40) 



where we used a slightly different definition of the current 
noise. This definition will be useful when examining the 
non-stationary case. Due to this definition the time de- 
pendance on t' in the symmetrized current noise is only 
present in the oscillator's quantum mechanical expecta- 
tion values. 

The general expression for the current noise in Keldysh 
space reads 

+ {x{t'))goi{t + t',t')+ 
+ {x{t + t'))gi^{t + t',t')+ 

+ iDit + t',t')gn{t + t',t')]. (41) 
where gij{t,t') is given in App. [A| 



1. General expression for the current noise 

In this section, we only make use of time translation 
invariance in the g^j (t, t') functions, the oscillator is taken 
as non-stationary. Details of the calculation can be found 
in App.jBj The final result we obtain for the symmetrized 
current noise to second order in the tunneling amplitudes 
reads 
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+ (i(t')> toh cos(77) 



{xit')) totiism{r/) 
toil 



a+{uj, V) + - {a+{Lj + n,V) + a+{u} ~ fl, V)) 
a- (cj, V)-]^ {a- {uj + n,V) + cr" (cj - VL, V)) 



{Pit')) 



ti 



in(r;) [a-{u + n,V) - cr"(w - n,V)] -icos{r]) [a+ {lo + n,V) - a+ (lu - n,V)] 



+ {x{t')x{t')) -j- \a+{Lo + n,V)+ a+{uj - n, V) 



2m 



{{Ht'),m}) ^ [^+(^ + n,v)- a+{uj - n, V)] } , 



(42) 



where we separated the real and imaginary part using the 
relative phase rj between the tunneling amplitudes and in 
addition introduced 



± 



2 



coth /? 



2 



(43) 



We want to note the important aspect of the current 
noise that it possibly can have a complex valued charac- 
ter which we discuss later. The gained expression is quite 
lengthly but provides us with the full quantum mechani- 
cal non-equilibrium characteristics of the current noise in 
the Markovian as well as in the non-Markovian regime. 
We made no assumptions on the state of the oscillator, 
which now gives us the possibility to identify momentum 
properties of the nanomechanical resonator using the cur- 
rent noise spectrum Ssym{^,t'). In the next section, we 
discuss this new possibility of a p-detector which involves 
measuring a complex valued current noise. 



2. Complex current noise and the p-detector m the 
non-stationary case 

The expression in Eq. ( |42[ ) allows for a comparison with 
results obtained in Ref. 1 161 where it was possible with a 
phase of 77 = 7r/2 to determine the momentum of the 
oscillator. In Ref. ilBj an Aharonov-Bohm setup allows 
the tuning of the relative phase rj. The full current noise 
spectrum there is proportional to the position spectrum 
Sx{i^) which is peaked aX uj = ±f2 in the case of 77 = and 
in the case of ry = 7r/2 is proportional to the momentum 
spectrum Sp{uj) showing peaks at w = ±f2. This peaked 
structure of the current noise spectrum is a fourth order 
effect, as we will see and discuss later when dealing with 
the fourth order corrections to the current noise. 



As one can see from Eq. ( 42 ) , already the second order 



current noise allows to determine the expectation value 
of the oscillator's momentum and in addition to that of 



the anticommutator even if the phase 77 = 0, i.e. 

the Aharonov-Bohm setup becomes obsolete in our case. 
The signature of the oscillator's momentum p in our case 
is however different than the one in Ref. [ini Instead of 
the peaked structure, we find a kink-like structure which 
stems from the fact that we deal with second order per- 
turbation theory. 

In order to understand how one can use this to identify 
the momentum, we have to understand the meaning of 
a complex current noise. As stated in Ref. [53 a complex 
valued current noise is in principle a measurable quantity. 
To have a relevant measurable quantity we would have 
to average the time dependent current noise Ssym{oJ,t') 
over the measurement time AT. We could do this in the 
following way 



sym V / 



1 

AT 



AT/ 2 



dt'Si%,iLO,t'). (44) 



-AT/2 



Since the time dependance of the current noise is only vis- 
ible in the expectation values of the oscillator's variables, 
it is important for the actual measurement to consider the 
time scales which are involved. If the measurement time 

AT is less than the time scale of the oscillator (l/fi), 

—(2) 

the measured time averaged current noise S'sym(i:Li) will 
be time-dependent. If however, the oscillator undergoes 
multiple oscillation cycles during the measurement time, 
the current noise will be time-independent. In this case 
we could as well take the oscillator to be in a stationary 
state. For a damped oscillator the times scales on which 
the damping happens have to be taken into account as 
mentioned already in Sec. |IVB[ 

We conclude with remarks on the interesting non- 
stationary case, where we can also take 77 = with- 
out losing the information on {pit')) and in addition 
obtain information on {{x{t'),p{t'y}) . We intend to 
give an idea of how to access the information on {p{t')) 
and {{x{t'),p{t')Y) available through the complex current 
noise. 

The expectation value of {p{t')) with respect to num- 
ber states or a linear combination of them, will al- 
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ways vanish when averaging over time according to 
Eq. (44). However, this is different for coherent states 
|a) = cxp(— |q;P/2) X^^o'^"/"'' !"■)' "^^isre we can write 
a — \a\exTp{id) with |a| being the amphtude and S the 
phase of the coherent state, respectively. The time aver- 
aged expectation values {p{t')) and {{x{t'),p{t')}) with 
respect to \a) yield 

(p(t')),„ = sin(17AT/2) , (45) 



and 



{{x{t'),P{t')})a. 



Ml 

nAT 



sm{2S) sin(17AT) . (46) 



For short measurement times AT < 1/fi we can write 

^hm^ (p(t'))a. - V2mn\a\ sm{S) , (47) 



2|apsin(2(5) 



(48) 



We sparate the current noise Ssym{co) ~ •^sym r(^) + 

— (2') 

•^sym li^) ^^^^ ^^^^ imaginary part where we observe 
that the imaginary part S^^'^ j(a;) only contains informa- 
tion on the oscillator's momentum {p(t')) and the anti- 
commutator {{x{t'),p{t')}) 



?(2) 
sym 

1 



_i(a.) = 27rp2e2{ 



toti\a\ sm{6) [a+{Lo + n,V)- a+{uj - n, V)] 

^ 3in(2(5) \a+(uj + n,V)- a+(u - fl, V)] ] . 
2mil J 



(49) 



The phase 6 of the coherent state now allows for a de- 
termination of the oscillator's momentum. For S = 7r/2, 
the signature in the imaginary part of the time averaged 

— (2) 

noise ' j(w) stems only from the oscillator's momen- 
tum. The signal in the non-Markovian regime is more 
pronounced than in the Markovian regime, see Eq. (49). 



3. Current noise in the stationary case 

Contrary to the non-stationary case we now also as- 
sume time translation invariance in the oscillator corre- 
lation function, i.e. D{t,t') = D{t — t'). One can see 
that the calculation in the stationary case goes along the 
same lines as in the non-stationary case. The only dif- 
ference will be that oscillator expectation values are now 
taken at time t' — 0, i.e. we encounter for instance (i(0)) 
instead of {x{t')). 

When interpreting the result for the stationary case we 
have to keep the constrains on oscillator expectation val- 
ues in mind. These constrains mentioned in Sec. IIV B II 



lead to vanishing expectation values of the anticommuta- 
tor {{x,p}) and vanishing expectation values for (i) and 
(p) . The current noise to second order is then equivalent 
to the ones previously obtained in Refs. [T7| and 1201 cf- 
Eq. (C.4) in Ref.i20iwith 72 = (x) = (p) = 0. 



C. Current noise to fourth order in the tunneling 
amplitudes 

We now turn to the investigation of the fourth order 
current noise. Since the fourth order perturbation theory 
involves a large amount of terms we use a diagramatic ap- 
proach. In the case of the fourth order current noise we 
restrict ourself to the stationary case for simplicity. An 
overview of all contributing terms in the non-stationary 
case is given in App. [C] In what follows we give a short 
explanation of the diagramatics. From Eq. (35) it be- 
comes obvious that S^"^^ contains fermionic expectation 
values which have the form 



■^ii,i2a3,i4('''l;'^2,T3,T4) = 

-f (3 ^ 4) - (2 o 4) + h.c.'^ . 



(50) 



The index ij determines whether we are dealing with /3o 
or ^1, ij G {0,1}. It is only necessary to evaluate the 



first expectation value in Eq. ( 50 ) using Wick's theorem. 



the other ones follow as indicated by (3 O 4), (2 o 4) 
and hermitian conjugation. This first term leads to 

X [Gi{ti,T2)Gi{t4,T3) - Gl{Tl,T3)Gl{Ti,T2)] X 

X [Gr{T2,Ti)GriT3,Ti) ~ GriT3,Ti)GriT2,T4^)] ■ (51) 

Similar to Ref. |T7| we use a diagrammatic representation 
for the expression in Eq. (50). As an example, we show 



the diagrams emerging from the expression in Eq. (51) 
in Fig. |3] and explain the components of the diagrams. 
Fermionic Green's functions of the reservoirs are repre- 
sented by solid lines and vertices are depicted by a dot, 
labeled with a time variable and a Keldysh index indicat- 
ing on which branch of the Keldysh contour the time lies. 
An integration over internal times ti and T2 is implicit. 
In addition we also have to sum over the two internal 
Keldysh indices k and I. We recognize that we have to 
deal with two types of diagrams: diagrams which con- 
sist of one closed fermion loop (diagrams in the lower 
panel of Fig. Isl) and diagrams consisting of two closed 
fermion loops/bubbles (diagrams in the upper panel of 
Fig. [3]) . These two different types of diagrams give very 
different contributions to the current noise which we will 
discuss below. We include the oscillator correlation func- 
tion D(t,t') in diagrams by a wiggly line connecting two 
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Gr(T3,T4) 



Grin, n) 



G,(T3,ri) 





Gi{ti.T3) 



Giin.Ti) 



GAn,T3) 



FIG. 3. Diagramatic representation for Eq. ( |51|l where we 
omitted the factor Pi-^fii^Pi^jS*^, note that Eq. (51 1 only con- 
tains fermionic Green's functions. 



vertices. In Fig. [4] we give an example of diagrams in fre- 
quency space containing one oscillator correlation func- 
tion. Here, integration over the two internal frequencies 
wi and u}2 as well as summation over the internal Keldysh 
indices k and I is implied. In frequency space the differ- 
ence between the closed loop diagrams, b) in Fig. [4] and 
the bubble diagrams, a) in Fig. |4] becomes clear: for the 
closed loop diagrams, the oscillator line always appears 
under integration of an internal frequency, whereas for 
the bubble diagrams, there is no integration over the os- 
cillator line. This is the reason for the different kind of 
contribution to the current noise of closed loop and bub- 
ble diagrams. As we will show below, bubble diagrams 
lead to peaks in the current noise, whereas closed loop 
diagrams lead to the afore mentioned kinks in the cur- 
rent noise. We reduce the number of diagrams by only 





FIG. 4. Examples of diagrams containing one oscillator cor- 
relation function D. In panel a) a bubble diagram is de- 
picted where the oscillator line is independent of an internal 
frequency. Panel b) shows a closed loop diagram where the 
oscillator line always appears under integration of an internal 
frequency. 

keeping contributions ~ igij and ~ tf which are the only 
finite coirtributions in the stationary case. Details are 
given in App.[C] This allows us to write the current noise 
S''''-'(t3, T4) for the further analysis as 



5(4) ( 



T3,T4) ^ S'j^\t3,T4) + S'}^\t3,T4)- 



(52) 



where S includes the bubble diagrams, S includes closed 
loop diagrams and D indicates the number of oscillator 
lines in the diagrams. The final result we obtain is valid 



for an arbitrary relative phase 1] which goes beyond the 
result obtained by Wabnig et al. in Ref. [171 This fact 
allows us to study the p-detector in fourth order pertur- 
bation theory. 



1. Results for S^^\tz,T4), S^^\tz,T4), 5'^D(T3,'r4) ana 



In the following, we sum up the different types of dia- 
grams, bubble type diagrams as well as closed loop dia- 
grams, both then can be integrated exactly. 

First we consider all diagrams containing only one os- 
cillator line, these diagrams are all proportional to t\ 
and depend on 77. We find for the symmetrized frequency 
dependent current noise 



cos(?7)^ 
- sin(77)^ 



4,e^V^ -AeVa-{w,V) 



Sx{oj) 



uj^ - 2ua+{oj,V) 



Sx{uj) 



+ i [D'^iuj) + D'^iw)] cos(77) sin(?7)x 
X [AeVa+{uj,V) -2ioa-{uj,V)] 



(53) 



which is one of the main results of this paper. 

In the case of the closed loop diagrams containing one 
oscillator line, it is also possible to sum up all diagrams 
and integrate them exactly. The expression for S^'^]^ ^ (w) 
is rather lengthly, therefore we do not present it here. 

We find that the current noise signature of S^yl^^ ^ (w) 

is of the kink-like structure similar to S^^^ j-, (a;) . In addi- 

(2) 

tion to the kinks at \uj\ — \ V±il\ coming from S'^y^ ^(a;), 

^symui^) gives rise to extra kinks at |w| — \V\ and 

\lo\ = However, these contributions are only minor 

(2) 

modifications to the current noise floor •S'gy^ ^(a;). Ex- 
periments as in Ref. il9i focus on the current noise near 
the resonance frequency a; « for which •SsymoC'^) 
the most important contribution. Therefore, the dis- 
cussion of our result will focus on the contributions 
stemming from S^yl^^ ^i^jj) ■ These contributions posses 
a peaked structure, since the oscillator correlation func- 
tions D^/'^{uj) and the spectrum Sx{uj) is peaked around 

The other contributions to the current noise stem from 
diagrams containing two oscillator lines D{t,t'). These 
diagrams are all proportional to f3if3l/3i/3l — t\ and 
therefore independent of the relative phase 77 between 
and ti. Moreover, these current noise contributions are 
small compared to the ones containing only one oscillator 
line since t\ <^ t^ti- We however are mainly interested in 
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the possibility to detect the oscillator's momentum which 
depends on rj, for this reasons and the fact that they are 

small compared to sj^'^]^ d('^) '^^ include them in 
our discussion, nevertheless we state our result which wc 
obtain after summing up the diagrams 



^(4)^7,,,^ _ 



^DD 



e 
'2n 



doji (fc/)| 
-D'^\u:,) D'^i-Lj,) g^ii^uji) + c.)} . (54) 



The last integration in Eq. ( 54 ) can be easily done since 



the oscillator correlation functions are peaked at ±17. 
Our result is then in accordance with the one obtained 
by Wabnig et al. in Ref. [T71 where it is has been shown 
that these contributions to the current noise are peaked 
at w = — 217, , 2f7 in contrast to the contributions arising 



from Eq. (|53|. Similar to s''^\^^), s'^^^i^) is of the kink- 



like structure and therefore only modifying the current 
noise floor. 

We now address the current noise stemming from 
S'gy^ j-,(aj) for arbitrary rj and also arbitrary system pa- 
rameters. 



2. Current noise in the Markovian and non-Markovian 
regime for arbitrary 77 

In order to compare our result for the momentum de- 
tector with Ref. dni we investigate sj.'^]^ -^{uj) near the 
resonance (w ~ 17). We find 



Ae^V^ cos (77) 
1 

— ^ sin(77)' 
- cos(77) sin(77) 



1 - 
1 



6^2 



m 472^^172 + (a;2- 172)2 
[4eT/cr+(17, V) - 217(7- (17, V)] 



Sp{uj) + 
(55) 



where cr^(^, V) is given in Eq. (431 and Sxii^j) and Sp{Ld) 
near resonance are given by 



S'x(i^) 



472^,172 + (c.2 _ 172)2 



(56) 



with X = x,p. The above expression is valid for the 
Markovian as well as for the non-Markovian regime. The 
relevant information about the oscillator can now be 
gained form the current noise spectrum. 



We take two different ways of evaluating the expec- 
tation values (a;2^ and (p^)- In the first one we use 
number-states which lead to expectation values {x"^},^ — 
{2n + l)/ml7 and (p2)^ — {2n + 1) ml7, where n denotes 
the oscillator's number of quanta. 

Since we also could imagine, as already explained in 



Sec. IV B 1 two equilibrium bathes, the tunnel junction 
and an external heat bath, to which the oscillator cou- 
ples, we can assign an effective temperature T^g to the 
oscillator which obeys 7tot T^s = 7o rcnv+7+ Tjunc, where 
7tot = 70+7+ is the total damping due to coupling to the 
junction (7+) and an external heat bath (70). The exter- 
nal heat bath is at temperature T^nv and the junction's 
temperature is given by Tjunc = eV/2kB- The oscillator's 
expectation values in this thermal regime are then give 
by (^2) = 2 /CB Toff /to 172 and (p2) ^2mkBT^G. 

For both cases, the thermal case and the number-state 
one, it is convenient scaling the current noise S'symoC'^) 
with e/o = 27rpQe2i;§tT+(17, V). We furthermore introduce 
dimensionless constant /i,/2,/3,/4 which are defined in 
the following way 



7tot 

7+ 
eV 



-^cnv 



17 

Ji 

7tot _ 
eV 



17 

hJ2 



(57) 

(58) 
(59) 
(60) 



/i can be interpreted as an overall quality- factor. The 
ratio 70/7+ = (/2 — 1) leads for G]1,2[ to a stronger 
coupling to the tunnel junction and for /2 > 2 to a 
stronger coupling to the external heat bath. The param- 
eter /a distinguishes the non-Markovian (/a e]0, 1]) from 
the Markovian regime (/a ^ 1). The last parameter /4 
quantifies the temperature Tdv of the external bath wit 
respect to the applied bias V . 

We now compare the signal of the position detec- 



tor S. 



(4) 

X— dot 



(a;) 



the momentum detector S, 



(4) 



0) to the signal of 



(4) 



'p-dotV"^/ — ^sym,D('^;'7 — 

7r/2) and later their dependencies on the parameters 
fi at resonance lo « 17. Assuming a high qual- 



ity resonator we take -^/l — 72^^/172 « 1 in Eq. (55 1. 
We call Qx = cr-(17, V^)/(4eFml7 (a;2^) quantum cor- 
rections to the x-detector current noise, arising from 
the non- vanishing commutator similarly we call 

Qp = (2toct"'"(17, y))/ (p2^ quantum corrections to the 
p-detector current noise. We then find 



S. 



(4) _.r2 [1 ~ Q=c] <j(4) 
dot ^ ^JS 



Qp] 



(61) 



and conclude that in the Markovian regime where /s 3> 1 
the signal of the position detector is always larger than 
the signal of the momentum detector. Whereas in the 
non-Markovian regime we have a stronger signature of 
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the noise part showing the momentum signature of the Figs. [6] and [7] show the current noise spectrum around 
oscillator. In the foUowing, we investigate in more detail a; w il in the Markovian and the non-Markovian regime 
the current noise of the x- and p-detector. for different couplings and environment temperatures. 



3. The x-detector 



10000 



From Eq. (551 one can see that for 77 = mod tt we 
recover the position detector result as in Refs.[T31[Ill and 
ITTI with peaks in the current noise spectrum at oj — ±51. 
Since we calculate the symmetrized current-current cor- 
relator, the current noise is symmetric in ui. The sign 
of the signal is given by the sign of [1 — Qx] which for 
an oscillator in the thermal regime depends on /2 , /s and 
/4, for an oscillator in number-state n it depends on n 
and /a only. We stick to a thermal resonator, noting 
that as in Ref. 16 for the p-detector, here the quantum 
corrections Qx ~ I//3 can become large (compared to 
1) in the non-Markovian regime where /a < 1, leading 
to a sign change. The parameter regimes for a nega- 
tive/positive peak in the current noise are depicted in 
Fig. [5] as blue/red regions. The change of sign in the 




1.10 



FIG. 6. Noise signal at a; = SI in the Markovian regime for 
different values of /2, /a, fi and /i = 200 which sliows a peak, 
due to the presence of the oscillator. 




FIG. 5. The current noise peak at lj — fl for the a;-detector 
as a function of /2, /s, /4, where we took /i — >■ 00. The blue 
region shows the parameter regime where the peak is negative, 
for parameter combinations lying in the red region, the peak 
is positive. 
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FIG. 7. a;-detector current noise signal at a; = f2 in the 
non-Markovian regime for different values of /2 , /a , A and 
/i = 200. The signal is weaker than in the Markovian regime 
(compare to Fig. [6|. In this regime however, it is possible to 
see a change of sign of the signal, depending on the parameters 
f2, fi- The inset illustrates this sign change. 



signal, depends on the environment temperature Tenv, 
the coupling to the environments /2 and heavily on the 
bias voltage and therefore /s. Deep in the Markovian 
regime the sign change is hard to achieve, only if /2 is 
very large and the bath temperature T^nv is very low, 
meaning that heating of the oscillator can be compen- 
sated by strongly coupling to a cold environment. In the 
non-Markovian regime the quantum corrections Qx can 
become large more easily and due to the lower signal in 
the non-Markovian regime for the x-detector, the sign 
change is more pronounced. 



4- The p-detector 

For the cases rj — -k /2 mod tt in Eq. ( [SS]), the result of 
the momentum detector as stated in Ref. 1 161 are extended 
to the non-Makovian regime. Due to the fact that the 
quantum corrections Qp are in the first place larger than 
Qa;, the peak in the current noise spectrum stemming 
from the oscillator has a negative sign for 77 = 7r/2. How- 
ever, it is also possible to change the sign by adjusting the 
parameters /2,/3,/4 on which Qp depends. In Fig. [8] we 
depict the regions with a negative sign blue and the ones 



12 



with a positive sign red. Changing the sign of the cur- 
rent noise signal in the p-detector case is easier to achieve 
over a wide range of parameters, as compared to the x- 
detector, even deep in the Markovian regime (/a ^ 1). 
Figure [9] shows a summary of the p-detector current noise 
in the Markovian and non-Markovian regime for different 
parameters /2,/4, respectively. 



lu 




FIG. 8. This figure shows the parameter regimes for the 
p-detector where the current noise peak at oj = S7 has a neg- 
ative sign (blue) and positive sign (red). Compared to the 
x-detector the sign change is possible for a wider range of 
parameters /2, /s, /i; we took /i — >■ 00. 



5. Detection of number states 

With the above made definitions we can map the oc- 
cupation number of the oscillator to experimentally ad- 
justable parameters, as for instance the environment tem- 
perature Tenv and the bias voltage V, similar to the ap- 
proach in Ref. [ITl This allows us in general to determine 
which state n the oscillator is in. This mapping is in- 
dependent of 77 and depends only on the dimensionless 
parameters fi in the following way 

2/21/31/4-21/31/4-/2 + 1/3! 
n = — . (62) 

The dependance of n on adjustable parameters is de- 
picted in Fig. [To] In the Markovian regime the oscillator 
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FIG. 9. The left panel shows the p-detector current noise 
in the Markovian regime for various parameters /2 , fi and 
/i = 200. The right panel shows the current noise in the 
non-Markovian regime for various parameters /2, /4. In both 
cases we can easily achieve a sign change in the signal. 




FIG. 10. Number of quanta n on the oscillator in depen- 
dence of experimentally adjustable parameters, for dominant 
coupling to the external heat bath (upper panel, with /2 = 5) 
and dominant coupling to the tunnel junction (lower panel, 
with /2 = 1.2). 

is only in its ground state for low environmental tem- 
peratures, since the applied bias voltage is heating the 
oscillator. In the non-Markovian regime, we can have a 
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higher environmental temperature for the oscillator being 
in its ground state. 

VI. CONCLUSION 

We have studied the finite frequency current noise of a 
tunnel junction coupled to a harmonic oscillator. In our 
work, we go beyond the Born approximation (because we 
calculate the noise up to fourth order in the tunneling am- 
plitude) and beyond the Markov approximation (because 
we do not restrict ourselves to the regime eV/Ml ^ 1). 
For a non-stationary oscillator, we have shown that the 
finite frequency current noise of the detector can be com- 
plex. This complex current noise can be used to obtain 
information about expectation values depending on x as 
well as expectation values depending on p. The former 
we call a;-detector signal and the latter p-detector signal. 

For the stationary oscillator, the finite frequency cur- 



rent noise is always real. Then, it is more complicated to 
get momentum information using a tunnel junction de- 
tector. An Aharonov-Bohm-loop setup is needed for this 
task. We analyze such a setup for the first time in the 
non-Markovian regime and thereby show how the x- and 
the p-signal can be used to determine the quantum state 
of the oscillator. 

Our analysis is an essential prerequisite to study how 
the quantum entanglement of NEMs2ll can be measured 
on the basis of tunnel junction detectors. This very in- 
teresting problem will be addressed in future work. 
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Appendix A: Details on the fermionic Green's functions 

Expectation values of the fermionic operators hi and ji can be expressed by the Keldysh Green's functions 
Gi,r[t,t') = -i(T,ci^r{t)clM')) as 



T,h,{t)j,{t')) = -ig^,{t,t') - -i [l3,p*Gr{t.t')Gi{t',t)-l3*(i,Gr{t\t)Gi{t,t')\ 
T,h{t)h,{t'))=g,j{t,t') = [l3,p*Gr{t,t')Gi{t\t) + l3*li,Grit',t)Gi{t,t')\ 

Tcm3Ai'))^QriM = [P^P■Gr{t,t')Gl{t\t) + j3*PJGr{t',t)Gl{tX)\ 



(Al) 
(A2) 
(A3) 



The Fourier transform of the function Qij{t^t') can be calculated in the usual way, yielding 



g-+{u)+ = 21, pI 



Z^^/')^— [-1 + coth(/3^— )] -f [1 + coth(/3^— )] 

[1 + coth(/3^— )] + [-1 + coth(/3^— )] 



(A4) 
(A5) 



duji 



\G^ (a;i-fa;)G; (wi) - (t^i)G, (wi -I- w)] 



(A6) 



where a (lo^V^ is given in Eq. (43 1. 
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Appendix B: Details of the second order current noise calculation 



The starting point for the calculation is Eq. (41), together with Eqs. (11 12 1 the symmetrized current noise in the 
non-stationary case can be written as 



dtt 



sm 



+ {xit')i{t')) cosint) [g^+it) + g+fit)] + ^'''^2^^^'^^ smint)g^+{t) + ^^^^'^Jf'^^ sin(m)g+f (O} ■ (bi) 



mil 



The further calculation is straightforward by using Eq. ( A4 1 and Eq. ( A5 1 . Finally, the current noise Ssyin{^-,t'^ in 
Eq. (Bl 1 can be written as 



{Pit')) [Gio^ (u^ + n)^ g^+ iuj~n) + g^V i^ + ^)- ^iV - ^)] + 



2m,n 
1 



+ - {x{t')xit')) [g^i+iu^ + + g^+{uj + g+-{uj + + g+-{u; - n)] - 



2mn 



{x{t')p{t')) [g^+{u + ^^)~ g^+{uj -n)]-^ W)m) [Qtii^ + ^)- Qtii^ - ^)] } 



(B2) 



The functions (T*(^, F), see Eq. (43), allow us to distinguishing the Markovian from the non-Markovian regime. For 
T we find 



'y-iLV) 



sgn(l^)e e|T/|>e 
sgn{V)e\V\ e\V\<£, 

'e\V\ e\V\>i 



where T here is the temperature of electrons in the leads. 



(B3) 
(B4) 



Appendix C: Details of the fourth order current noise calculation 



We first give the whole expression for the current noise to fourth order in the tunneling amplitudes containing the 
7\4-functions of Eq. ( 50 ) and oscillator operators 



5(4)(r3,T4) = -y^dridr2 { 



A^O,0,0,o(ti,T2,T3,T4) + 

+ A^0,0.0,l(n,T2,r3,T4) (x(t4)) + A^O,0,l,o(n,r2,T3,T4) {x{tz)) + 
+ A^O,l,0,o(Tl,T2,r3,T4) (x(t2)) + A4l,0,0,0(n,T2,T3,T4) {x{ti)) + 

+ A4o,0.1,l(n,T2,r3,T4) {TcX{t3)x{t4)) + A^O, 1,0,1 (^1 , T2 , r3, T4) {TcX{T2)x{Ti)) + 
+ A^O,l,l,o(n,T2,r3,T4) {TcX{T2)x{T'i)) + A^l,0,0,l(n,T2,r3,T4) {Tcx{ti)x{t4)) + 

+ A^i.o,i,o(Ti,T-2,r3,T4) {Tcx{ti)x{t3)) + A^i,i,o,o(ti,t-2,t-3,T4) {Tcx{ti)x{t2)) + 

+ A4o.i,i,i(n,T2,r3,T4) {Tcx{T2)x{T'i)x{Ti)) + 1,0,1,1 (ti , T2 , T3 , r4) (T^ x(ri)i (t3)£ (T4)) + 

+ -Ml.l,0,l(n,T2,r3,T4) {TcX{Ti)x{T2)x{Ti)) + 1, 1,1,0 (ti , T2 , T3 , r4) {TcX{Ti)x{T2)x{T:i)) + 

+ Xi,i,i,i(ri,r2,r3,T4) (T, i(Ti)x(r2)x(r3)£(T4)> } + /" dn dT2 (t, Htun{Ti)i{r3)'j (T,7?tu„(r2)/(T4)) . (CI) 
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A first reduction of terms in Eq. (CI I is done by only focusing on the stationary case. This allows us to drop 
terms which are proportional to {x{t)) of Eq. (CI) and only keep terms proportional to i') a,nd D{t,t') D{t" ,t"'). 
Unlinked diagrams which appear in this expression are canceled by the term which is always of the bubble type. 
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